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Abstract 

We address the question of the existence of magnetic domains in a recent scenario 
for the creation of primordial magnetic fields. We write down the magnetohy- 
drodynamic equation in an expanding universe in the infinite conductivity limit 
and discuss its solution. No domains are found. 
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We have recently presented a scenario for the creation of a primordial magnetic 
field at the GUT scale fl|]. This field is due to nonperturbative quantum effects in the 
early universe. For a typical realistic GUT such as the supersymmetric SU(5) model 
we found a field of the order B ~ 5 x lCU 8 /i 2 where // is the renormalization point 
with /i ~ T GUT ~ 10 15 GeV. 

Once the field is imprinted on the charged plasma it will remain because the early 
universe is a very efficient conductor with conductivity inversely proportional to the 
collision cross section. For details we refer to the Appendix of Turner and Widrow ||, 
where the conductivity has been estimated. It then follows that whatever magnetic 
flux there existed in the early universe, it is frozen in the plasma. Diffusion begins 
relatively late and starts from small scales. 

In order to estimate the magnetic field today one needs to take into account the 
fact that in the universe there exists causal horizons. Naively one could argue that 
each horizon "bubble" (i.e. a bubble with a radius equal to the causal horizon at the 
GUT time H~ x ~ 10~ 28 cm) constitutes a magnetic domain. Thus one should average 
over all random domains, which now have a size of about 1 m, in order to obtain 
the average magnetic field today. If this were so, then the value of B at the time of 
protogalaxy collapse would very likely to be too small to be of interest for the observed 
galactic magnetic fields. 

In the present note we shall, however, show that the picture of a causal horizon 
bubble as a magnetic domain with randomly pointed field, such as in a ferromagnet, 
is not correct. Instead the proper physical situation is that when two bubbles meet, 
a new "joint" bubble is formed with the magnetic field lines frozen to the plasma 
particles. After a short time the joint plasma looses all memory of its origin, and so 
does the magnetic field, too. The situation is analogous to having two containers filled 
with identical gases which are separated by a wall. Once the wall is removed the joint 
system has no memory of the original two separate containers. 

When the electric conductivity is very large the only electric field that can exist 
is induced by the motion of charged particles in the magnetic field. Here it is very 
likely essential that the field B is very small relative to the scale /z 2 ~ T 2 , as is the 
case in our ferromagnetic universe model [|]]. If the field is large, say of the order of 
T 2 , the mean free path of the charged plasma particles is large in comparision with 
their radius of curvature in the magnetic field. In such a case conductivity cannot 
be estimated from the drift of the charged particles and in general it depends on the 
magnetic field itself. 

We use the same conventions as in ref. 0. In particular, we adopt the conformal 
metric 

dr 2 = a(r]) 2 (-dr] 2 + dx 2 + dy 2 + dz 2 ) , (1) 
where r] is the conformal time. The advantage of the conformal metric is that electric 



and magnetic fields a treated on an equal footing with 
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a 2 B 3 (and cyclic), F i = —a 2 Ei, (etc.) (2) 



The standard cosmological model described by Eq. ([I]) is now assumed to be valid 
on the average. We treat the metric as a background field not affected by the presence 
of the magnetic field itself. In the relation physical length scale = ax (comoving length 
scale) the comoving length scale related to a charged plasma particle is assumed to 
be fixed (as in the standard cosmological scenario) only on the average. Similarly, 
the charged particles are only assumed to be at rest at a given Friedman- Robert son- 
Walker time in an average sense. Thus, the particles are associated with a velocity 
field v(x, t). With the definitions (|j) one of the Maxwell equations reads 

da 2 B „ - 

— = -V x a 2 E , (3) 

where the spatial derivative is with respect to the comoving coordinates. Now, because 
of the very large conductivity E must be an induced field, so that 

E = -v x B . (4) 

Hence Eq. (pj) becomes 

<9a 2 B , 

— = V x (v x a 2 B) , (5) 

which is the magnetohydrodynamics equation for infinite conductivity (see || for dis- 
cussion of this subject) in an expanding universe. 
Eq. (H) can be written as 

^^ = a 2 (B-V)v-(v-V)a 2 B-a 2 BV-v. (6) 

In order to proceed we need to know the divergence of v. In the nonrelativistic limit 
(ignoring a) it is given by the usual continuity equation 

^ + V-(nv) = 0. (7) 
In general relativity Eq. (0) is replaced by 

("^-^ + ^v^) = 0, (8) 

where n is the density of charged particles which is conserved, and U^ 1 is the four- 
velocity. In the conformal metric we obtain 

- V • v = — L- |-(a 4 *7 n) + -^v • VnU° . (9) 



Inserting this in Eq. (H) we obtain 



d ( a 2 B \ ( d \ a 2 B 

' — + v-V 



drj \a A U°n J \drj J a A U°n 

(a 2 B- V 



a A U°n 



(10) 



in an adiabatically expanding universe n ~ a~ 3 . To find how U° scales, let us note 
that in the FRW metric f/^ w = dt/dr = 1, so that in the conformal metric U° = 
dn/dr = l/a(rj). Of course, in our case this result is valid in the average only, and 
v = U /U° is locally non- vanishing (although (v) = 0). Thus we find that 

a A nU° ~ const. (11) 

Eq. flUf ) can be solved. Let X(t/i) be the position of a plasma element at the time 
771, and let x = x(X, 772) be the position of the same plasma element at the time 772. 
Then Eq. (|10|) is solved by 

a(^) 2 Bi(x,r/ 2 ) _ a(r/i) 2 B j (x,7/i) dxi 



a(772) 4 n(r/ 2 )f/ (772) a^n^U* dX 3 

Because of Eq. ( PD we see that a 2 B changes by the factor dxi/dXj. 

Another way of expressing the solution (|T2| ) is the following. Let 5xi be an infinites- 
imal displacement vector carried by the fluid (in comoving coordinates). If initially 
Sxi(0) is parallel to a 2 (0)5,(0), 

6xi(0) = 00(0)2.0,(0) , (13) 

with c a constant, then the solution (0) implies that the two vectors Sx and a 2 B move 
together with the plasma, and at later times 

5x0?) = ca(r]) 2 B(ri) (14) 

In the standard cosmology the comoving displacement vectors satisfy 

(Sxi(0)) = (Sxiir,)) , (15) 

and hence a(w) 2 Bi(r]) is a constant. This means that the magnetic flux measured in 
the physical coordinates is constant, which is the well known result. 

Let us consider B 2 , using the special solution ( p^ -flU). It is clear that 

a{nYB 2 {n) = a^B 2 , 77) = a(0) 4 B 2 (x(0), 0) (16) 



statistically. Eq. flllf) is identical to what happens for e.g. the particle density, for 
which one finds 

a(77) 3 77(x(77), 77) = a(0) 3 7i,(x(0), 0) . (17) 



Now, returning back to the original problem, namely whether different GUT hori- 
zon bubbles form domains or not, we observe that a 4 B 2 behaves like the particle 
density a 3 n. Of course, the density a 3 n does not have any domains. If two horizon 
bubbles collide, a new joint bubble is simply formed, and the information about the 
original bubbles is lost (after all, in equilibrium entropy is maximized). Exactly the 
same thing happens with respect to a 4 B 2 . If in the initial set-up there are magnetic 
domains, these are rapidly washed out by the motion of the plasma (gas). 



So far we have used the special solution (pf)- (|T4]) . However, we may also take the 



general solution (|12|)> which by virtue of Eq. (|TT| ) can be written as 

- 

a( % ) 2 B i (x( % ), % ) = a(r /l ) 2 B J (x(7 ?1 ),r ?1 )^ . (18) 
From a statistical point of view it is natural to assume isotropy so that 

<j|>=V (19) 

In other words, if we have an original displacement SXj at r\ — T]\, then at a later time 
r)2 the corresponding displacement 5xi can in principle go in any direction. However, 
the original j-direction is the only preferred direction, so that if j turns into any i ^ j 
then it is equally likely that the resulting displacement is positive or negative. Hence 
the statistical sum is zero. Eq. flT9| ) is of course also valid when we consider the 
standard cosmology as a statistical model. 

Thus Eqs. ([Tj|) and ([19]) again lead us to the conclusion that a 4 B 2 is (statistically) 
constant. They also give rise to the strong result that 

a(r ?2 ) 2 (B,(x(r /2 ),r ?2 )) = a(7/ 1 ) 2 (B i (x(7/ 1 )^i)} , (20) 

which demonstrates that, apart from a well understood redshift factor, the field B is 
statistically constant. Eq. (EDf) shows that the evolution of a 2 B is exactly the same 



as the evolution of an invariant density such as a 3 n (see Eq. (|I7])) having no domains. 
The same result holds for a 4 (B 2 ). It follows that the magnetic energy Eb scales like 
l/a{rf) and hence Ebci remains a constant during the course of the evolution of the 
universe. Consequently the magnetic energy always corresponds to a magnetic field 
of 10 _8 T 2 . Physically this is so because in the infinite conductivity limit there is no 
dissipation. 

We have thus made it plausible that in the framework of the standard cosmology 
considered in a statistical sense there are no magnetic domains. Therefore one may 
estimate the magnitude of the primordial magnetic field today as was done in our 
previous paper ffl with the result that 



5L ay - (10" 14 G) 2 



(21) 



This certainly has the right order of magnitude to act as the seed field for the galactic 
dynamo mechanism. 

Let us end by reiterating that our argument very much depends on the fact that in 
our case the primordial magnetic field is relatively weak. If the 5-field is strong, of the 
order of 0(T 2 ), conductivity would depend on B and simple magnetohydrodynamics 
would no longer be valid. In the construction of Vachaspati |4]], relying on the fluctu- 
ations of the Higgs fields at the electroweak phase transition, the field is locally 10 8 
times larger than in our present model. Hence it is quite possible that such a strong 
magnetic field does not freeze in the plasma and domains may exist in that case. 
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